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Abstract. Limb-darkening is a fundamental constraint for modeling eclipsing binary and plan- 
etary transit light curves. As observations, for example from Kepler, CoRot, and Most, be- 
O |e come more precise then a greater understanding of limb-darkening is necessary. However, limb- 
D ' darkening is typically modeled as simple parameterizations fit to plane-parallel model stellar at- 
^0 , mospheres that ignores stellar atmospheric extension. In this work, I compute linear, quadratic 
I ■ and four-parameter limb-darkening laws from grids of plane-parallel and spherically-symmetric 
' model stellar atmospheres in a temperature and gravity range representing stars evolving on the 
Red Giant branch. The limb-darkening relations for each geometry are compared and are found 
I I to fit plane-parallel models much better than the spherically-symmetric models. Assuming that 
limb-darkening from spherically-symmetry model atmospheres are more physically representa- 
^0 five of actual stellar limb-darkening than plane-parallel models then these limb-darkening laws 
will not fit the limb of a stellar disk leading to errors in a hght curve fit. This error will increase 
C>) with a star's atmospheric extension. 

^ , Keywords, stars: atmospheres, stars: fundamental parameters, (stars:) supergiants 
1. Introduction 

^— I Stellar limb- darkening is the observed change of intensity from the center of the stellar 

^ disk to the observable edge, where the intensity decrease is due to the geometric pro- 
jection of the line-of-sight relative to the radius of the star. This effect is an important 
challenge for the interpretation of observations of binary stars (e.g. Claret 2008), and 
■^llj^ ■ planetary transits (e.g. Knutson ct al. 2007, Croll et al. 2011), as well as interferometric 
' (e.g. Wittkowski et al. 2004) and microlensing (e.g. Zub et al. 2011) observations. Typ- 
. ically, limb-darkening is treated as a parameterization or relation as a function of the 
' cosine of the angle formed by the radius and line-of-sight, called ji to simplify analysis. 
T-H , Stellar atmosphere models and binary/transit observations are complementary tools for 
^ understanding limb-darkening and stellar astrophysics in general because observed limb- 
• ^ i darkening can help constrain models. There are a plethora of articles describing different 
/\ ' limb- darkening relations (Al-Naimiy 1978, Wade & Ruchinski 1985, Claret et al. 1995, 
^ Claret 2000), limb-darkening coefficients from predicted intensity profiles for a number 
" " of stellar atmosphere codes such as Atlas, and Phoenix (e.g. Howarth 2011, Sing 2010, 
Claret & Hauschildt 2003, and different fitting methods (Wade 1985, Heyrovsky 2003, 
2007, Claret 2008). In this work, I focus on a small number of limb-darkening laws and 
compare predicted fits for intensity profiles from plane-parallel and spherically symmetric 
model stellar atmospheres. In the next section, I describe the stellar atmosphere code 
and three limb-darkening laws of interest: a linear, quadratic, and four-parameter (Claret 
2000). In Sect. 3, I present results of the fitting the limb-darkening laws using model 
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atmospheres, and how the errors of the fit depend on assumed geometry. I summarize 
this work in Sect. 4. 



2. Stellar atmosphere code and Limb-darkening Laws 

For this analysis, I use a new Fortran 90 version of the Kurucz Atlas code (Lester & 
Neilson 2008). The code computes opacities using opacity distribution functions, atmo- 
spheres are assumed to be in local thermodynamic equilibrium, and hydrostatic equilib- 
rium. Each atmosphere model outputs intensity profiles as a function of wavelength, for 
an equal spacing of /i for 1000 points. Typical calculations for the Atlas code compute 
intensity profiles for 10 - 17 /i-points. The program computes models for either plane- 
parallel or spherically-symmetric geometries, where the plane-parallel model is described 
by two fundamental parameters such as T^ff and log 5, while spherical models require an 
additional parameter such as stellar mass. Neilson &; Lester (2008) fit model intensity 
profiles to interferometric observations from Wittkowski et al. (2004) and predicted sim- 
ilar fundamental parameters as those authors. Also, Neilson & Lester (2011) predicted 
limb- darkening coefficients for a specific limb-darkening law from spherical models and 
compared to results for microlensing observations from Fields et al. (2003) and found 
better agreement than the authors did using plane-parallel models and spherical models 
that had intensity profiles clipped to remove the extended limb. 

I have computed approximately 2000 model stellar atmospheres in spherical symmetry 
for the parameter range Tcs = 3000 - 8000 K in steps of 100 K, log 5 = — 1 - -I- 3 in steps 
of 0.25 in cgs units, and M — 2.5 - 10 Af© in steps of 2.5 M©. Plane-parallel models are 
computed for the same values of Toff and log g. 

For this work, I compute least-squared fits to three laws: 
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1 — a(l — ji) Linear, (2-1) 

1 - 6(1 - ^i) - c(l - /i)2 Quadratic, (2.2) 
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- ^ d,(l - /i'/2) Four Parameter, (2.3) 

i=l 

where intensities are computed in the Kepler white light passband. All fits are computed 
using least-square fitting of the limb-darkening coefficients. The quality of the fit may be 
measured in a number of ways; here, I test the quality of the fit of limb-darkening laws 
by checking how well they conserve stellar flux, AF/F ~ (FModci — -FLaw)/-FModci- 



3. Results &: Summary 

In Fig. [TJ I show the predicted intensity profiles for a Toff — 4000 K and logy = 2 
model atmosphere for both geometries. There is a significant difference between the 
model intensity profiles such that the intensity near the limb of the spherically-symmetric 
model atmosphere is much smaller than the plane-parallel model atmosphere. The plane- 
parallel model does not appear to go to zero, though the equation of transfer suggests 
that the intensity in the limit as /i — > then I{fi) (Mihalas 1978). This may be 
an issue with the resolution of fi. The plane-parallel model is also much better fit by 
the limb-darkening laws than the spherically-symmetric model atmosphere because the 
spherically-symmetric intensity profile is more complex. 

The relative difference between the model intensity profiles and limb-darkening laws 
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Figure 1. Intensity profiles for a plane-parallel (left) and spherically-symmetric (right) model 
atmosphere with T^s = 4000 K and \ogg — 2, along with the best-fit limb-darkening laws. 

predicted stellar fluxes are shown in Fig. [51 The fits to plane-parallel model atmospheres 
appear to be better, the average error is < 5% for the linear law, and is much smaller 
for the other laws. The quality of the fit for plane-parallel model atmospheres is also ap- 
parently independent of effective temperature. The results for the spherically-symmetric 
model atmospheres are strikingly different. The flux errors are much larger, 5-10% for the 
linear law, 0-20% for the quadratic law, and 0-5% for the four parameter limb-darkening 
law. The difference in fits due to model atmosphere geometry suggests a significant prob- 
lem for understanding stellar limb-darkening. It is reasonable to assume that a spherical 
geometry is a more physical representation of an actual star than a plane-parallel stel- 
lar atmosphere, hence people should be hesitant for using limb-darkening coefficients 
generated from plane-parallel model atmospheres. It also suggests that these are not 
ideal limb-darkening laws to use and it may be necessary to develop new limb-darkening 
relations for the future. 
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Figure 2. Relative difference of between model and fit stellar fluxes for the three limb-darken- 
ing law in the Kepler passband, linear (upper- left), quadratic (upper- right), and Claret (2000) 
four-parameter (bottom) laws for plane-parallel and spherically-symmetric model stellar atmo- 
spheres. 
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Discussion 

R.E. Wilson: Your starting explanation of why the intensity goes to zero at the limbs 
for a plane parallel case is not correct. Actually the intensity does not go to zero at the 
limb. The mistake came from neglect of emission along the line of sight it is not just 
an attenuation problem, but has both a source factor and an attenuation factor in the 
intensity integral. If one looks into an infinite uniform region, the received intensity is 
not zero, but is the intensity characteristic of the regions temperature. 

I. HuBANY: To Bob Wilson's comment: The intensity does not indeed have to go to 
zero at the limb, but such a case is not covered by a 1-D plane-parallel treatment of 
the transfer equation anyway because in this case the medium is infinite with no natural 
boundary condition. Comment to the talk: The Eddington factor (the ratio of the K- 
moment and the mean intensity) is not necessarily equal to 1/3 in the plane-parallel 
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case. Such quantity is usually called a variable Eddington factor, and depends on depth 
and frequency; it goes to 1/3 only deep in the atmosphere. 

A. Prsa: I understand why one would use an analytic model for a Mandel-Agol type 
approach, but perhaps the systematic error from the simple fit may be avoided simply 
by linearly interpolating (or looking up) /(/x). 

Ph. Stee: I did not understand why you used the first lobe of the visibility function to 
fit the LD instead of the second lobe, especially since you may fit the first lobe with a 
uniform disk? 



